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ABSTRACT 

ON ' Recent results on standard solar models are reviewed. I shall summarize briefly 



On 



three of the themes that I stressed at the Neutrino '92 Conference: 1) Different solar model 



codes give the same answers when the same input data are used; 2) Improved calculations 

o ■ ....... 

of standard solar models include helium diffusion, the Livermore radiative opacity, the 
' meteoritic iron abundance, and a variety of other corrections; and 3) There are a few 

basic rules that should be followed in using standard solar models. At the Neutrino '92 
Conference, I reviewed in more detail the recent work on standard solar models by Marc 
Pinsonneault and myself. This work has by now appeared in print (Rev. Mod. Phys. 64, 

OO ■ 

Os| 1 885, 1992, hereafter Paper I, and ApJ Letters, 69, 717, 1992, Paper II). Therefore, there 



is no need for me to repeat the details here. 
1. DIFFERENT CODES GIVE THE SAME RESULTS FOR THE SAME INPUT DATA 
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a: 

Since the first accurate solar model was derived in 1962, a number of different codes have been developed 
to calculate precise solar neutrino fluxes It was pointed out on page 78 of Neutrino Astrophysics, Cambridge 
University Press (1989) that the results of nine different codes developed between 1962 and 1989 were 

•rH , 

consistent to within about 10% for the 8 B neutrino flux when the same input parameters were used or when 
corrections were made for the differences in input parameters. References to the various calculations are 
given on the just-cited page 78. 

Some authors have not been aware of this previously-established convergence of calculated fluxes and 
attributed differences between fluxes calculated using some contemporary solar models to different stellar 
physics. Historically, a few examples of program errors have been discovered, but these have all been 
relatively minor in their effects. In fact, the contemporary differences are essentially all due to different 
choices in input data. This result was demonstrated in a step-by-step comparison of recently-derived models 
in Section 4 of Paper I. 

The final result of this detailed comparison is that, with the same input parameters, the standard solar 
models of Bahcall and Ulrich; of Sienkiewicz, Bahcall, and Paczyhski; of Turck-Chieze, Cahen, Casse, and 
Doom; of Sackmann, Boothroyd, and Fowler; and the current Yale model (Bahcall-Pinsonncault) all predict 
event rates for the chlorine experiment that are the same within 0.2 SNU. The Sackmann, Boothroyd, and 



Fowler results are also apparently in good agreement with the other calculations but Boothroyd has informed 
me of approximations that were embodied in their stellar evolution code that prevent a precise quantitative 
comparison. 

2. IMPROVEMENTS IN SOLAR MODELS 

A number of improvements have been made in the input data for standard solar models over the past 
four years. I will first discuss the effects of helium diffusion, radiative opacity, and iron abundance, all of 
which have comparable effects (about 1 SNU) on the predicted event rate in the chlorine experiment. 

The new physics that is most difficult to include in stellar evolution codes is the diffusion of helium 
with respect to hydrogen in the solar interior. Convenient formulae were derived for this effect by Bahcall 
and Loeb, but these formulae still cause complications in the stellar evolution codes because diffusion mixes 
spatial and time derivatives. So far, the only precise solar models including this effect were presented in 
Paper I. Helium diffusion increases the predicted event rate in the chlorine experiment by about 0.8 SNU. 
That is to say, if you want to compare two otherwise similar solar models in their predicted event rates, 
you must add SNUs to the predictions of the model that does not include helium diffusion: 0.8 SNU for the 
chlorine experiment and 5 SNU for the gallium experiment. 

A major improvement in the radiative opacities was achieved through the calculations of Iglcsias and 
Rogers, whose results are incorporated in the now widely used Livermore Opacities. These opacities give 
better agreement in comparing stellar model predictions with observations of Cephei pulsations, RR Lyrae 
masses, double mode Cepheids, double mode 8 Scutti stars, Li abundance in the Hyades cluster, the depth 
of calculated convective cores, and the luminosity-effective temperature relation for O and B stars. The use 
of Livermore rather than Los Alamos opacities increases the calculated event rate by about I SNU for the 
chlorine experiment. 

There has long been a major uncertainty in the predicted event rates due to the different values inferred 
observationally for the abundance of iron on the solar surface. In Paper I, we adopted-following recent 
reanalyscs of the reliability of the photospheric abundances-the lower iron abundance obtained by studying 
meteorites rather than the higher abundance obtained by analyzing the solar photosphere. This choice 
reduces the event rate by 1.3 SNU. 

Paper I also describes a large number of less important corrections including more accurate nuclear 
cross sections for several fusion reactions, a more detailed electron screening calculation (although a small 
correction may still need to be made in this area, see the discussion by Turck-Chieze in this volume), and a 
few programming improvements. 

The predicted rate for the chlorine experiment given in Paper I is 

< fa >ci, theory = (8.0 ± 1.0) SNU, (1) 
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where I have reluctantly quoted here-in deference to my experimental friends-an approximate la error 
(instead of my preferred "total theoretical error" , which includes 3a errors for all measured quantities) . The 
difference between the experimental rate of 2 SNU presented at this conference and the theoretical rate given 
in Eq. (1) constituted for two decades the "solar neutrino problem." 
The predicted rate for the gallium experiment is 

^(<M Ga - 131.51J SNU, (2) 

again approximate la errors. The GALLEX result is about 2a less than the calculated rate; the SAGE 
result is more than 2a less than the theoretical value. 
The calculated flux of 8 B neutrinos is 

<H 8 B) = 5.7(1 ± 0.15) x 10 6 CHi^V 1 , (3) 
which is approximately twice the measured Kamiokande rate reported at this conference. 

3. HOW TO USE STANDARD SOLAR MODELS 

Physicists are used to comparing the results of different experiments that measure the same or related 
quantities. We arc not used to comparing different solar models. However, for the next few years we 
must continue to use the predictions of accurate solar models to interpret what we have learned from solar 
neutrino experiments. With the next generation of experiments (SNO, Super-Kamiokande, ICARUS), precise 
measurements of the 8 B energy spectrum and the neutral current to charged current ratio will allow us to 
decide if there is new physics required independent of solar model calculations. But, these results will 
presumably not become available until sometime in 1996 or later. 

What should one do to be an informed consumer of solar models? Here is a list of questions that you 
should ask about solar models. 

(1.) Does the model yield the same results (to about ±0.1 SNU for the chlorine experiment) 
as other standard models with the same input data? 

Paper I verifies that this condition is satisfied for all of the contemporary models discussed at this 
conference. 

(2.) Have the authors used the best published values (and the quoted experimental errors) on 
all measured quantities? 

Most people calculating stellar models do not bother to continually upgrade their calculations as more 
accurate measurements of, e.g., nuclear quantities become available. Some models discussed at this conference 
were not adjusted for the recent thorough reanalysis by the Cal Tech group of the 8 B production cross section 
(with a derived error). An exportable nuclear energy generation and neutrino production subroutine that 
makes this easy to accomplish is can be obtained from the author. 
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(3.) Does the model include the best-available theoretical input data? 

Of particular importance are radiative opacities, element diffusion, and equation of state. Of the models 
discussed in this conference, only the models in Paper I include the improvements due to helium diffusion 
and to the Livermore radiative opacities. Most of the models include in the adopted equation of state the 
dominant correction caused by the Debye-Huckel effect. 

(4.) Are the errors well defined and equal to the published errors determined by the experi- 
mentalists? 

For anyone interested in the subject of errors on solar models, Section 8 of Paper I contains a current 
discussion of the best-estimates of the errors and a comparison of how different errors have been estimated. 
Some authors have quoted subjective uncertainties that are larger than the values given here and in Paper 
I. I believe that the well-defined errors derived in Paper I and in this report are conservatively large. The 
reader is invited to consider Section 8 of Paper I and to judge on the basis of the detailed discussion whether 
the quoted uncertainties are indeed conservative. 

4. CLOSING COMMENTS 

We need more solar neutrino experiments. It is possible that we as a community may have been lucky 
and have stumbled onto some important microscopic physics while trying to test accurately stellar evolution 
theory. But, the available results depend upon a small number of experiments with redundancy only in 
the gallium experiments and with a continued reliance on the standard solar model calculations. Unless 
the gallium results converge together on a value that does not exceed 80 SNU (cf. chapter 11 of Neutrino 
Astrophysics) even with the experimental errors, we will probably have to wait for the next generation of 
experiments (SNO, Super-Kamiokandc, Borcxino, and ICARUS) before we can prove that we have learned 
new physics or new astrophysics. 
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